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Abstract 

The upside-down simple harmonic oscillator system is studied in 
the contexts of quantum mechanics and classical statistical mechanics. 
It is shown that in order to study in a simple manner the creation 
and decay of a physical system by ways of Gamow vectors we must 
formulate the theory in a time- asymmetric fashion, namely using two 
different rigged Hilbert spaces to describe states evolving towards the 
past and the future. The spaces defined in the contexts of quantum 
and classical statistical mechanics are shown to be directly related by 
the Wigner function. 
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1 Introduction 

In this work we will study the motion of a particle subject to an upside-down 
simple harmonic oscillator potential 

V{q) = ~muj\^ (1.1) 

in the contexts of quantum mechanics and classical statistical mechanics. 
The aim of the paper is threefold: 

1. -We solve the problem using a nonconventional technique. In fact, we 
find that the cvohition of the system can be entirely expressed in terms 
of idealized states that decay or grow exponentially (usually called Gamow 
vectors). In order to give mathematical meaning to these states we are 
forced to work in the framework of a rigged Hilbert space (RHS). In quantum 
mechanics the Gamow vectors are then defined as generalized eigenvectors of 
the Hamiltonian operator with complex eigenvalues, being these the poles of 
the scattering matrix when it is extended to the complex plane. The Gamow 
vectors can be used in 'generalized spectral expansions' similar to those found 
by the Gel'fand-Maurin theorem, but in order to use them we have to define 
two different test function spaces, that we shall call and A similar 
result is found in classical statistical mechanics. 

2. - We will prove that with this mathematical structure we have set the 
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basis to introduce a time asymmetry in the theory. Let us explain this fur- 
ther. The vectors in can be expanded as hnear combinations of the 
decaying Gamow vectors while the vectors in can be expanded as linear 
combinations of the growing Gamow vectors. Since every physical state (ei- 
ther quantum or statistical mechanical) must decay both towards the future 
and the past, it is mathematically sound to represent the state of a physical 
system when it evolves towards the future (i.e. from an initial condition) by 
a vector 0+ G and the state a system when evolving towards the past 
(i.e. going to a final condition) by a vector G When considering the 
scattering of particles by the potential barrier, this division between initial 
states and final states is easily explained, but we think that it is useful for 
studying other phenomena, like the decay towards the equilibrium position. 
Since both spaces and are dense in the corresponding Hilbert space 
Ti,, at this stage there are no empirical results that could help us decide if 
the conventional representation of the state of a physical system by a vector 
in Ti is better than the representation of the same system by a vector in $_|_ 
or Both representations give the same empirical results, so both are, in 
a sense, correct. The difference is at a mathematical level (in the topologies 
used, etc.) and so we can use the representation we find more suitable to the 
description of the physical facts. 

Now, by selecting two different mathematical structures to represent a 
physical system when evolving towards the future and towards the past we 
have introduced the basis for a time asymmetry in the mathematical de- 
scription of time evolution. In fact, there are only two causes for asymmetry 
in nature: either the laws of nature are asymmetric or the solutions of the 
equations of the theory are asymmetric. E. g.: the laws governing the weak 
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interaction are asymmetric while the solutions of the theory are asymmetric 
in the case of spontaneous symmetry breaking. 

Time-asymmetry is not an exception. Thus, if we want to retain the time- 
symmetric laws of nature the only reason to explain the time-asymmetry of 
the Universe and its subsystems is to postulate that the space of solutions 
is not time-symmetric, namely to use the second cause for asymmetry. So 
the proper way to solve the problem is simply to define a realistic time- 
asymmetric space of admissible physical solutions $+; namely to restrict the 
space of initial conditions. If K is the time inversion operator, this space will 
be time-asymmetric whenever K : i— namely time inversion changes 
the physically admissible solutions in $4. into a space of inadmissible solu- 
tions $_ different to the previous one. If Ti is the usual space of solutions of 
the theory (Hilbert space in quantum mechanics or Liouville space in clas- 
sical statistical mechanics) and we select properly, then mathematically 
speaking we have introduced a Gel'fand triplet C C 

A Reichenbach branched system is perhaps the most realistic model for 
an irreversible Universe, i. e. a set of irreversible processes such that each 
one begins in an unstable state produced by another member of the system 
and it eventually ends in an equilibrium state This set of processes, all 
of them beginning in a non-equilibrium state, defines a global arrow of time 
in the Universe. The problem is that the whole branch system must begin 
in a global unstable initial state which has no explanation. This unstable 
initial state would be the initial cosmological state of the Universe. It is 
qualitatively shown in ref. |0] that the expansion of the Universe can be 
the agency that produces this initial unstable state. Using the method of 
this paper, we have found the same explanation but in a quantitative way 
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[0], endowing the Universe with a (global) space of admissible solutions $+. 
Therefore, we think that the time-asymmetry is not given by the system 
itself, namely our upside-down oscillator; instead it must be connected with 
the global arrow of time of the Universe. In fact, our system is really a 
member of the Reichenbach branch system in such a way that the unstable 
initial condition of our oscillator is necessarily produced by another member 
of the branch system. Thus, as an initial condition is given, it is admissible 
only if it belongs to a (particular) space $+. Since the Rigged Hilbert space 
and the conventional Hilbert Space formulations are both equally correct, we 
can select the space better qualified to define time-asymmetry. 

3.- As this procedure has been applied in the past to some unstable 
quantum mechanical systems (simple potential scattering problems 0] and 
Friedrichs' model 0) and chaotic classical statistical systems (Renyi maps 
[P] and Baker's transformation 0), in this paper we extend this technique 
to the simplest unstable system: the upside-down harmonic oscillator. Fur- 
thermore, we show that the two sets of spaces defined in quantum mechanics 
and classical statistical mechanics are connected in a simple manner. This is 
the first result of this kind that we know of. 

The organization of the paper is as follows. In sec. 2 we give the principal 
properties of the rigged Hilbert space formulation of quantum mechanics. 
We define the regular state space and the generalized state space and we 
define generalized eigenfunctions. In sec. 3 we work out the motion of a 
particle subject to the potential ( |1.1| ) in classical mechanics, introducing 
the canonical variables we shall use throughout the paper. In sec. 4 we 
study the same problem in the context of quantum mechanics. We define 
the growing and decaying Gamow vectors and using these objects, we also 
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define the two RHSs that represent states evolving to the future and states 
evolving to the past. In sec. 5 we undertake the same task in the context 
of classical statistical mechanics. We find the "statistical Gamow vectors" , 
namely generalized density functions in state space that represent idealized 
growing and decaying states. Once again, using these objects we define the 
two RHSs that represent states evolving towards the past and towards the 
future. In sec. 6 we show that the structures defined in sections 4 and 5 are 
connected by the Wigner function. Finally, in sec. 7 we draw our conclusions. 
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2 The rigged Hilbert space formulation of quan- 
tum mechanics 



In the traditional (von Neumann's) formulation of quantum mechanics, a 
physical state is represented by a vector in a Hilbert space H. and physical 
magnitudes (observables) by linear selfadjoint operators acting in it. We will 
deal in this paper with the unidimensional motion of a particle in a potential 
field, so the Hilbert space we should work in is isomorphic to L^(IR). More 
precisely, in the position (|g)) representation the particle's state is represented 
by a normalized wavefunction ilj{q) G whose modulus squared gives the 

probability density of finding the particle in the position q. Observables arc 
then represented (in this representation) by selfadjoint operators in L^(m). 

This formulation contains certain idealizations, since once we interpret 
the wavefunction as a probability amplitude, the only physical requirement 
is that it must be square integrable. But a Hilbert space is a complete topo- 
logical space, with respect to a particular topology, namely the one obtained 
from its scalar product. The assumption that every vector in this Hilbert 
space represents a physically realizable state cannot, by no means, be justi- 
fied by empirical facts, since a topology (infinite limits, continuity) has no 
physical meaning. It is just a mathematical idealization with which a theo- 
retical physicist works in order to formulate a theory. So we can take another 
mathematical idealization and formulate a theory in a different mathematical 
environment. For instance, we can take another topological space, which is 
itself a subspace of the Hilbert space, and associate with the vectors in such 
a space the states of the physical system. 

This is what the rigged Hilbert space (RHS) formulation of quantum 
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mechanics does 0]. To define a RHS we take a topological vector space (en- 
dowed with a nuclear topology) $ with a continuous scalar product defined 
in it. As we shall see, when the Hilbert space is the space of square inte- 
grable functions, the test function space $ is chosen as the set of functions 
of fast decrease or Schwarz functions S = 5(k) or some subspace of it. By- 
completing the vector space with the topology given by the scalar product 
we get a Hilbert space Ti such that $ C 7^. If we consider the set of all 
antilinear functionals, continuous with respect to the nuclear topology, we 
get another vector space called the dual space of $ and denoted $ ^ . We will 
denote the value of the functional F G $^ on the vector G $ by (</)|F) and 
its complex conjugate by {F\(f)). Due to Riesz's lemma, the dual space of 7i 
is H itself, so we get (using the fact that the nuclear topology is stronger 
than the Hilbert space topology) the Gel'fand triplet 

<^>c7^c<^>^ (2.1) 

We associate the vectors in the space $ with the physical states of the 
system in consideration, thus it is usually called the regular state space. 
The observables are associated with continuous (essentially) selfadjoint linear 
operators in $. 

This formulation has some advantages over the conventional one. The 
first advantge is that in the traditional formulation of quantum mechanics 
a vector is a class of Lebesgue square integrable functions differing in a set 
of measure zero while in the RHS formulation a vector is, usually, just one 
continuous and infinitely differentiable function. 

The second advantage is that every observable is well defined, since they 
are continuous in $. This means in particular, that even when the operator 
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is unbounded in Ti. it is well behaved in $, which is contained in the domains 
of all operators of interest. In the case we will be studying, these include the 
position, momentum hamiltonian operators. 

The third advantage is that every essentially selfadjoint continuous linear 
operator in a RHS has a complete set of generalized eigenvectors in $^ 
(the generalized state space) with their eigenvalues in the spectrum of the 
operator, a result proved by Gel'fand and Maurin |[. 

By a generalized eigenvector of an operator A in a RHS we mean a func- 
tional I Fa) G $^ such that 

(AVIFa) = A(</.|F,) y<j)E^. (2.2) 

Since this is a generalization of the definition of an eigenvector in finite 
dimensional vector spaces, the number A is called a generalized eigenvalue 
corresponding to the eigenvector \Fx). What the Gel'fand-Maurin theorem 
states is that given an operator, the set of generalized eigenvectors spans 
the regular vector space $. More precisely, if A is the spectrum of the 
aforementioned operator A, then the scalar product between two vectors 
0, G $ can be expressed as 

(0,V^)= / df^iX){^\F,){F,\i;) (2.3) 

J A 

where /i is a certain integration measure. We see then that in the RHS 
formulation of quantum mechanics, Dirac's notation is fully justified. Even 
though the generalized eigenvectors do not belong to the Hilbert space, they 
are defined as antilinear functionals in In the same way, the opera- 

tor itself can be expressed as a linear combination of the same generalized 
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eigenvectors, namely 



A = [ dfi{\)\\\){\\. (2.4) 

J A 

There is one more advantage that this formulation has over the traditional 
one: we can define generalized eigenvectors of an essentially selfadjoint op- 
erator with eigenvalues that do not belong to the (Hilbert space) spectrum 
of the operator. In the general case, such a spectrum is a closed subset of 
the real line but, as we shall show below, we can find in some cases eigenvec- 
tors with complex eigenvalues. This will allow us to define Gamow vectors, 
namely generalized eigenvectors of the hamiltonian operator with non zero 
imaginary eigenvalues. As we shall see, this fact implies that the evolution 
of these vectors is exponential, either growing or decaying, depending on the 
sign of the imaginary part of the eigenvalue. 

These new generalized eigenvectors are useful to study the temporal evo- 
lution of the regular vectors representing physical states, if we define the 
RHS so that we can find new expansions like ( |2.3| ) containing them. We 
will find that in order to do so, we must define two different RHSs, that we 
shall call and The first one will correspond to the representation of 
physical systems when they evolve towards the future, since its vectors will 
be expressed in terms of the decaying Gamow vectors, and the second one 
will correspond to the representation of physical systems when they evolve 
towards the past, since they wil be expressed in terms of the growing (or 
decaying towards the past) Gamow vectors. 

To finish this section we will mention the different rigged Hilbert spaces 
with which we will work in this paper . The first one is the one constructed 
from the space of Schwarz class functions S, composed of all infinitely differ- 
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entiable functions of a real variable such that together with all their deriva- 
tives vanish at infinity faster than the inverse of any polynomial. The RHS 
obtained is 5 C L^(m) G . In this case, the functionals in are called 
tempered distributions. This is the RHS most theoretical physicists work in, 
since in it the position and momentum operators are continuous 0]. 

The other two RHSs we shall use are constructed from subspaces of S, 
so that this last property is maintained. These subspaces are the space /C 
of infinitely differentiable functions of a real variable with compact support 
and the space Z of Fourier transforms of functions in /C, that is isomorphic 
to the space of entire functions of fast decrease. The RHSs we obtain are 
then /C C L2(m) c /C^ and Z C L\r) C Z^. 
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3 The classical upside-down oscillator 



The system we are going to study is one of the simplest unstable ones, namely 
the motion of a particle in the presence of a potential of the form (|1.1|) . The 
hamiltonian function of the system is H = — ^muj'^q^. 

As in the study of the harmonic oscillator, it is convenient, at this point, 
to adimensionalise the dynamical variables. In order to do that, we must 
take the natural scales of length, momentum and time defined through the 
physical parameters available: m, % and uo (we include % since we will deal 
in the next section with the quantum case). These scales are, respectively 



h/muj , a/ mujh and l/cj, so we are making the transformations 



-^q , p ; — (3.1) 



and 



H^^H , t^ujt. (3.2) 
riuj 

The relation between the Hamiltonian and the adimensional q and p vari- 
ables is 

H=\{p'-q'). (3.3) 

To solve the equations of motion, it will be helpful throughout the paper 
to work in another couple of canonical variables 



v = ^{p + q) , u=^{p-q) (3.4) 
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obtained through the generating function F(q,v) — —\v'^ + \/2qv — \q^. 
Expressed in these new variables, the hamiltonian reads 



H — vu 



(3.5) 



and thus, the equations of motion are 

dv dH . du dh . , 

dt du ^ dt dv ' 

These equations are uncoupled and they have as solutions for any initial 
condition vq, Uq 



In phase space, the trajectories of the particles are hyperbolic and the v 
and u axes arc the corresponding asymptotes. The point f = 0, ?i = (or 
= 0, p = 0) is a point of unstable equilibrium. If Hq — vqUq ^ then the 
particle decays (gets away from the barrier region) both towards the past 
and the future. The directions parallel to the v axis will be called "dilating 
fibers" while those parallel to the u axis will be called "contracting fibers" . 



U — UqC 



(3.7) 
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4 The quantum upside-down oscillator 

4.1 The 1^;) and \u) representations 

As we have seen in the previous section, the evolution of the system is ex- 
pressed in a simple way when we use the variables v and u. In quantum 
mechanics, a canonical transformation is associated with a change in rep- 
resentation. Instead of using the representations |g) and \p) of generalized 
eigenfunctions of Q and P, we will use the representations \v) and \u) of 
generalized eigenfunctions of the operators V and U, defined through 

V = 1={P + Q) f/=-L(P-g). (4.1) 

Since these operators are the quantum representations of canonically con- 
jugate variables, they satisfy the commutation relation 

[V,U]=iL (4.2) 

The spectrum of these operators is the whole real line |p , and the trans- 
formations from the |g) representation to the \v) and \u) representations are 
given by 

{q\v) = (2vr2)-i/4e^(v^'"'-^'/2-V2) (4,3) 

{q\u) = (_27r2)-i/V(v^"5+^'/2+«V2)^ ^4_4) 
Due to (|4.2|), we get the relation 
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4.2 Eigenfunctions of the Hamiltonian 
4.2.1 Real eigenvalues 

The Hamiltonian of the system is 



H = \P^ - \q' = \{VU + UV) (4.6) 



namely, the symmetric version of ( |3.5| ). The eigenvalue equation for the 
eigenfunctions of this Hamiltonian with eigenvalue e in the \v) and \v) rep- 
resentations reads 



v^0e(^) = (^e - ]^)(t)e{v) , u^^,{u) = (-ze - ^)^.(u). (4.7) 
Formally, the solutions of these equations are 

(j),{v) = av''-'/^ , Mu) = Pu-''-'/^ (4.8) 

but care must be taken, since these expressions are only defined for pos- 
itive values of v and u respectively (in fact, as we have seen in sec. 2, they 
are not functions but distributions). 

Actually, there are two linearly independent solutions of the equations 
(^4.7|) for each value of e, due to the degeneracy of the Hamiltonian. These 



independent solutions can be chosen as ITD 
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{v\t + {v)) = ^e{vy^-'i^ , {v\e-{v)) = ^e{-v)\vr-'i^ (4.9) 

or 



{u\e+{u)) = ^e{u)u--''' , {u\e-{u)) = ^9i-u)\ur^~'/^ 
y In y 



(4.10) 

The generalized functionals ( [4.9[ ) represent the ideahzed scattering out- 
states, representing particles leaving to the left and the right respectively, 
while the generalized eigenfunctions ( fl.lO] ) represent the scattering in-states 



representing particles entering the scattering region from the left and the 
right respectively |Tl|, |13|. Both of these sets of solutions form complete and 
orthonormal sets in iS, in the sense discussed in sec. 2. 

By calculating the scalar product between the in-states and the out-states, 
we obtain the scattering matrix, whose coefficients are in this case of the form 



ra, 13 



5^.(e)=^.(e)r(^-ze) /i, z/ = +, - (4.11) 
where the f^v{e) are entire functions of e. 

4.2.2 Eigenfunctions with complex eigenvalues 

As can be seen from the above formula, the scattering matrix when extended 
to the complex plane has an infinite number of imaginary poles located at 
Zn = —i{n + 1/2) where n is a nonnegative integer. The Gamow vectors will 



16 



have as generalized eigenvalues these numbers or their complex conjugates. 
Instead of taking the conventional way to obtain the expressions for these 
Gamow vectors (namely by analytical extension of the scalar product to the 
complex plane we will take a shortcut and find them in an heuristic 

way. 

We will consider the solutions of the eigenvalue equations ( [4.8| ) when 
taking e = 2;„ in the first one and e = z^ in the second one. Then we get the 
functionals 



{v\n) = t;" , {u\h) = ^^^m" (4.12) 

v27m! 



where the multiplicative factors have been selected so that ([4.14]) below 
applies. Transforming with ( |4.5| ) we get 



(—1 

(u\n) = v^r5(")(M) , (v\n) = ^-^^("Vt;). (4.13) 

n! 

It is clear that these functionals are tempered distributions. By direct 
calculation, we can demonstrate the biorthonormality of the sets {\n)} and 
(Ifi)}, namely that 

{n'\n) = {n\n) = 6n,n'- (4-14) 

We will show now that \n) is indeed a generalized eigenfunction of H with 
complex eigenvalue Zn, namely that {H(j)\n) = z„(0|ra) Vvr G S. 
We have 



+ 00 



d(f){v) 



{H4>\n)= / dv[t{v^ + -<P{v))y 
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and integrating by parts in the right hand side, we get 

f+oo I I 

{H(j)\n) = -i dv{ (j){v) {n + (j){v)} = -i{n + -){(f)\n). 

In a similar fashion, we can demonstrate that \n) is a generalized eigen- 
function of H with eigenvalue ^ = i{n + 1/2). It is 

j-oo dv 2 n\ 

and usmg ■^{v^) = v^-^ + we get {H(l)\n) = 2„(0|n). 
Since \n) and \n) are generalized eigenvectors of H, their temporal evo- 
lution is easily calculated. In fact, we get 

e-^^» = e-("+i/2)*|n) , e~'^'\n) = e^^+^/^^'ln) (4.15) 

showing that the Gamow vectors are, indeed, vectors that would rep- 
resent idealized states that decay or grow in a perfectly exponential way. 
These functionals are more pathological (are "less physical" ) than the eigen- 
functions of the Hamiltonian with real eigenvalues ( [4.9|) or ([4.10|) ; they are 
clearly distributions and not regular functions and thus, cannot represent, 
by themselves, physical states. It will be shown, though, that they are useful 
in studying the temporal evolution of the regular states. 

To complete the presentation of these functionals, let us study their ex- 
pression in the |g) representation. In order to do that, we use ([4.4|) and 
get 
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where is just a numerical factor. Now, using the formula Hn{z 



9" a2+2A2^|^^^ where Hn{z) is the n-th Hermite polynomial, we find that 



{q\n) = a^e'^'/'H^{e-'^/'q). (4.16) 
In an analogous way, using (^73|) we get 



{q\n) = d^e-''^''/'H^{e'-/'q). (4.17) 
Restoring the variables' dimensions, we get 



{q\n) = C„e— «^/2^ff„(^^g) , {q\n) = de-^'"-^^/^^ff„(^^g). 
If we compare these generalized states with the eigenstates of the har- 



monic oscillator with equal frequency and mass [|I2|, we can see that the for- 
mer can be obtained from the latter by means of the transformation u t— > —iu 
in the case of \n) and the transformation (—> zcj in the case of |n). But 
these transformations applied to the potential of the harmonic oscillator turn 
it into the potential ( |1 . 1] ) and the same happens with the eigenvalues: they 
are transformed from discrete real eigenvalues to discrete imaginary ones. 
Once again we see (now in the |g) representation) that these generalized 
eigenfunctions cannot represent physical states; when calculating the norm 
squared of these functions we obtain that they diverge in the limit |g| i-^ oo 
as q . 
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4.3 Generalized expansions and the regular function 
spaces and 



We study now the use of the generahzed functions ( ^.12| ) in generahzed expan 



sions. Our first step towards this end is taking such that (f |0) = (p{v) G S. 
Since this wavefunction is infinitely different iable, we can define its Taylor 
expansion around v = 

= E ^^v- = T.{^\n){n\<j>). (4.18) 

n=0 n=0 

The equality in this formula is restricted to the values of f in a interval 
on the real line, namely inside the radius of convergence of the series, so the 
utility of this expansion is rather limited. Only when considered in a scalar 
product will it proof really useful. 

Our second step is to define the test function spaces 

$+ = e n/{v\^) eZ} = {^e n/{u\^) g JC} (4.19) 
= e n/{v\^) elC} = {tlje n/{u\^) e Z} (4.20) 

where K. and Z are the spaces introduced in sec. 2. We can then construct 
two different RHSs with these spaces, namely 

$± C 7^ C (4.21) 
Let us take (f) G $+; in this case the expansion (|4.18|) is valid for all values 



of V. Considering that \n) is a generalized eigenfunction of H with eigenvalue 
Zn and the definition ( [4.12| i), the time evolution of this vector is given by 
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(P{v,t) = (t;|e-'^V) = J2e-^''+^/^^\v\n){n\(f)) = e-'/^(j){ve-\0). (4.22) 

n=0 

Since $+ is a dense subspace in H, we get that ([4.22|) must be valid for 
any vector in Ti; this can be verified by direct substitution in Schroedinger's 
equation. If we consider the Taylor expansion of a vector in $_ in the \u) 
representation, we get that in that representation the temporal evolution is 
given by ^{u,t) = e*/^'?/'(ue*, 0). 

These two results can be seen as the fact that the wavefunction does not 
change its form with time; it just suffers a change in scale. If the wavepacket is 
initially concentrated in some value vo (resp. uq) then it will be concentrated 
at time t in foe* (resp. mqc"*). Namely, the center of the wavepacket follows 
the classical trajectory found in sec. 3. This is a consequence of Ehrenfest's 
theorem, since the potential is quadratic in position. Quantum effects come 
from the broadening of the wavefunction in the \v) representation and the 
narrowing of it in the \u) representation. 

Now, let us consider the scalar product between a vector 0_ G and a 
vector 'ip^ e $+, {(j)-,ip+) = J^P^ dv (j)_{v)il)+{v). Since 0-(f) G /C, it is an 
infinitely differentiable function of compact support; the integration limits 
can then be replaced by —a, a for some a. On the other hand, since ip+{v) is 
an entire function, the radius of convergence of its Taylor expansion around 
f = is infinite, and so we have 



(0_,^+)=/ dv<p.{v)Y.{h\ij+)vr (4.23) 

n=0 

The convergence of the series in the interval [—a, a] is uniform and, since 
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the function 0_ (v) is bounded in that interval, we can interchange the order 
of the summation and integration. Then we get 

oo 

(0_,^+) = ^(0_|n)(n|^+). (4.24) 

n=0 

We get a second expansion of this kind by taking the complex conjugate 
of this expression: 

oo 

(^+,</._) = E(^+l^)(^l</>-)- (4.25) 

n=0 

This shows that a vector 0+ G can be expanded (when acting as 
a functional in as = J2 and a vector 0_ G can be 

expanded as |0„) = 

Let us now turn to the physical meaning of these RHSs. As we have 
seen, if 0+ G the expansion ( [4.18|) is valid for all values of v. This means 
that we can represent a state as an infinite sum of decaying states. Then, if 
we want to study the decay of a physical system we can represent its initial 
condition by a vector 0+ G If on the other hand, we want to study the 
creation of a physical system then by representing the final condition with 
a vector ■?/;_ G we can think of this process as the growth of a linear 
combination of the exponentially growing states \n). 

In an experiment, we control the initial condition, wich evolves towards 
the future, and then compare it with a final condition by means of a mea- 
surement; thus, the quantities of interest are of the form (■?/'_, 0+ (t) ) = 
('0_ |e~'^*|0+), wich can be expanded as in ( [4.24| ) 



22 



(^_|e-^^*|0+) = ^e-("+i/2)*(</.^|n)(n|^+). (4.26) 

n=0 

It is evident from this formula that we can think that the initial con- 
dition is fixed and the final condition evolves to the past. The series in 
the formula converges always (for all values of t) as we have seen, but it 
only has physical meaning for positive values of t. For long times, only 
the first term of the series makes a contribution and then we can take 
{^p-, (f)+{t)) = e~*/^(0_|O)(O|^/'+) showing that in this regime the decay is 
effectively exponential, unless one of the coefficients (0-|O) or {0\ip+) van- 
ishes. We recover thus the lifetime already known from reference |TB|]. We 



remind the reader that equations ([4.24|) and ( ^4.26|) are exact, there are no 



approximations in the series. In our system, unlike the one considered in 
there is no background term; all the details of the evolution are found in the 
series expansions. 

Finally, let us consider time reversal in our formulation. From the inter- 
pretation we have given to the spaces and it seems that the time 
reversal operator should transform one of the spaces into the other, since it 
changes initial conditions into final conditions. This we shall show now. 

Wigner's time reversal operator is defined in the |g) representation as the 
conjugation operator, namely 



K : 10) 10') where 0'(g) = 0(g). (4.27) 
Let us take G $+. Then we have 
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ol/4 p+oo 



(4.28) 



21/4 



dqe 



-i{V2vq-g'^/2~v^/2) 



V27ri "'-oo 

By comparing ( [4.28| ) and ( [4. 291 ), we get 



0(g). (4.29) 



and considering (|4.19| ) and ( [4.20|) , 0' G , or 



(4.30) 



K : $1 I — > $_ 



(4.31) 



The inverse relation can be proved in a similar way. 
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5 The upside-down oscillator in classical sta- 
tistical mechanics 

5.1 Generalized eigenf unctions in classical statistical 
mechanics 

In classical statistical mechanics, a state of a physical system is represented 
by a density function p(t>, u) that gives the probability density to find it in 
a given region in phase space. This density function must satisfy certain 
requirements to have a physical meaning. First of all, it must be a positive 
function, since probabilities must be positive. Secondly, it must be integrable 
over the entire phase space. Finally, we ask the further condition that it be 
square integrable, since we want to calculate mean values of functions that 
will be themselves square integrable functions. In this case, we work in a 
Hilbert space. Just like we did in quantum mechanics, we will formulate 
the classical statistical theory in a rigged Hilbert space rather than in a 
Hilbert space, thereby introducing a time asymmetry in the representation 
of physical states. This method has been applied before to simple chaotic 
problems 0, 0. 

Let us consider the equation that gives the temporal evolution of the 
density functions, namely the Liouville equation 




(5.1) 



where the Liouvillian operator is defined by 



Lp = i{H,p}p,B. = ii 



OH dp OH dp 



(5.2) 



dv du du dv 
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( |5.1| ) is formally equivalent to the Scliroedinger equation, so the solutions 
to this equation can be found in a similar way. The Liouvillian is an (essen- 
tially) selfadjoint operator in the Hilbert space of square integrable functions 
of two real variables L^(]R^). We will find the generalized eigenfunctions of 
this operator that play a similar role as ( [4.13|) did in sec. 4. 

The Liouvillian operator is, in our case, Lp = i{u^ — v^) Taking v as 

the eigenvalue of L and proposing as a solution the product of a function of 
f by a function of -u (p = V{v)U{u)), we get 

uU-^U' - vV-^V = -iv. (5.3) 

It can be seen then, that both terms on the left hand side of this equation 
must be constants, that we shall call m and n, respectively: 

uU-^U' = m, vV~W = n. (5.4) 

We see that the equations for U and V are the same, so they will have sim- 
ilar solutions. A solution obtained by direct integration is V = Av^, U = 
Bu"^ where A and B are arbitrary constants. We get then that p = Cv'^u^ 
is a generalized eigenfunction of the Liouvillian operator with eigenvalue 
u = i{m — n) . As we did in sec. 4, let see what happens when we take m 
and n to be nonnegative integers; in this case the eigenvalue is imaginary. 
We will denote these eigenfunctions as 

\m,n) = — — -f"M™ ,h' = i{m — n). (5.5) 

They are merely polynomials in v and u. Once again, we find that this 
functions have no direct physical meaning; in fact, they do not satisfy any of 
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the properties we asked for the physical density functions. 
Using the relation 

we can see that V = 6'-^'\v) is another solution to ( |5.4[ 2) if we take 
n = — (n' + 1). We find then the following generalized eigenfunctions of the 
Liouvillian operator with their corresponding eigenvalues 

\m, n) = , u = -i{m - n) (5.6) 

|m,n) = ^ — ,u = i(m + n + l) (5.7) 
ml 

\m,n) = ^ — ;^(5(™)(M)f" ,u = -i{m + n+l). (5.8) 
nl 

These generalized eigenfunctions are biorthonormal by pairs, namely 

{m,n\rh'h') = {m,h\m' ,n') = Sm,m'Sn,n'- (5.9) 

All this functional are tempered distributions, acting on the space = 
5(r2). 

5.2 Rigged Hilbert Spaces and time asymmetry in Clas- 
sical Statistical Mechanics 

With the generalized eigenfunctions just found we can construct four different 
function spaces. Nevertheless, to find time asymmetry we must work with 
( |5.7|) and (|5.8|). It has to be that way, since we must treat in a different 



27 



way the contracting and dilating fibers, which in this case are the u and v 
variables. 

On the other hand, let us consider the temporal evolution of these gen- 
eralized eigenf unctions of the Liouvillian. We get 

|m,n(t)) = e(™+"+i)V,^) , |m,n(t)) = e-("^+"+^)*|m,n) (5.10) 

namely, the first ones correspond to idealized states that grow exponen- 
tially towards the future while the second ones represent idealized states that 
decay exponentially towards the future. These are the "statistical Gamow 
vectors" that we wanted to find. In a completely analogous way to what we 
did in sec. 4, using ( p.lO| ) we find that the solution to the Liouville equation 
is p(f , u, t) = p(t>e~*, ue*, 0). 

We define, as in sec. 4, the space of vectors that represent physical states 
when studying their evolution to the past and to the future, by 

^+ = {p{v, u) e L^/p e Z{v) ® )C{u)} (5.11) 

= {p(v, u) G LVp G }C{v) (g) Z{u)}. (5.12) 
Let us now consider the scalar product in L^(IR^) 

/+00 r+oo 
dv / du p{v,u)p' {v,u) 
-oo J—oo 

where the conjugation has no effect since the functions are real. Given 
p+ G ^'+, we get 
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where the coefficients of the series are infinitely differentiable functions 
of u, that vanish outside a given interval [—a, a] of the real line. Similarly, if 
p_ G \1'_ then 

P-(V,U) = r— (f , 0) 

where now the coefficients of the series are infinitely differentiable func- 
tions of V that vanish outside a different interval [—6, b] of the real line. 
Working as we did in sec. 4, we find 



m,n=0 



(P-,P+)= ^P-\'n^^n){m,h\p+). (5.13) 

m,7i=0 

In a restricted sense (when used in an expression between a function 
p+ G \1'+ and a function p_ G we get 



1= ^ |m,n)(m,n|. (5-14) 

TO,n=0 



Considering the complex conjugate of ( |5.13| ) we get 



1= ^ |m, n)(m, n|. (5.15) 

m,n=0 
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Just as in the quantum mechanical case, the justification for the dis- 
tinction between the mathematical representations of a physical state of the 
particle, when studying its future evolution and its past evolution, comes 
from the fact that in the former case we want the expansion in terms of de- 
caying states (I5.IOI 2) to be valid and in the latter case we want the expansion 



in terms of the growing states ( |5.10| i) to be valid. 



Let us now face the problem of time reversal in classical statistical me- 
chanics. Once again we shall find that the time reversed version of a function 
in \E'_ belongs in \E'+ and vice versa. In this case it is easier to see, since time 
reversal in classical mechanics is just the transformation 

q ^ q , p ^ —p , t —t 
and then the transformation for the variables f ad -u are 

V I— > —u , u —V (5.16) 

and so, aside for sign change, one of the variables is transformed into the 
other. Looking at the definitions of the spaces \E'-|- we get that \—>- $zp. 
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6 Connection between the classical and quan- 
tum cases 



To complete our study on the introduction of a time asymmetry in an upside- 
down simple harmonic oscillator system, we will show how the structures 
defined in quantum mechanics (sec. 4) and in classical statistical mechanics 
(sec. 5) are connected. This connection is given by the Wigner function. 
This function is defined from the density matrix p representing a quantum 
system by the formula 

F,(g,p) = (27r)-i j^^ dy{q - y /2\p\q + y /2)e'^y . (6.1) 

It gives an idea of the probability density in classical phase space if the 
system is represented by p [1^. In our case, we are considering pure states. 



so the density matrixes are p = \ip){'>p\ and then (|6.1| ) reads 



+00 



F^(g,p) = (27r)-i / rfy ^(g - y/2)^(g + y/2)e^P^. (6.2) 



Since we have been working mostly in the \v) and \u) representation, 
it will be useful to express the Wigner function in this variables. It can 



be shown [O] that the Wigner function is invariant under linear canonical 



transformations, then we can take 

/+00 
dyi,{2v-y)^{y)e'^^^y-^^\ (6.3) 
-oo 

Let us consider that G $+ and then '4){v) G /C. If the support of ip{v) 
is [—a, a], then the integration will be performed in this interval. We can see 
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that Wigner function is infinitely derivable in its two variables, since we can 
make the derivation inside the integration sign. 

Now, if we fix the u variable and look at the dependence on the v variable, 
then we can see that it vanishes if |f | > a. In fact, given that \2v — y\ > \2v\ — 
\y\ > a G [—a, a], there is no interval in the real line that contributes 
with a non vanishing term to the integral. 

On the other hand, let us fix the v variable. We see then that the Wigner 
function as a function of u is the Fourier transform of an infinitely differen- 
tiable function of bounded support, and thus, is a function in Z. We found 
then the relation 

IV') e $+ I — > F^{v, u) G IC{v) ® Z{u) = (6.4) 

In a similar way, we can get the relation connecting the two spaces that 
describe the physical systems when they evolve to the past 

G I — > F^{v, u) G Z{v) ® /C(m) = (6.5) 

This result is not unexpected, since in our case the potential is quadratic 
in the position. Under this circumstances, the equation that rules the tem- 
poral evolution of the Wigner function coincides with the Liouville equation 
|15[1 , so there must be a close relationship between the structures found in the 
classical statistical case and the density function we get from the quantum 
case. 
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7 Conclusions 



Wc found that in order to give meaning to Gamow vectors, two different 
rigged Hilbert spaces must be defined: one to represent states evolving to- 
wards the future from initial conditions and one to represent states evolving 
towards the past from final conditions. This distinction between initial con- 
ditions and final conditions must be thought of as an implementation of the 
global arrow of time of the Universe and not as a manifestation of an intrinsic 
irreversibility pertaining to the system. Nevertheless, the instability is the 
key to the implementation of the time asymmetry by causing the existence 
of Gamow vectors. 

By using this time-asymmetric formulation of the problem, we found that 
new "generalized expansions" can be used to make calculations in an easy 
way. For example, we found the solution of the Schroedinger and Liouville 
equations by means of the Gamow vectors. The fact that we restrict the 
space of vectors representing physical states has no empirical consequences 
since the test function spaces we work in are dense in the corresponding 
Hilbert space the conventional theory is formulated in. 

Even though our model has some characteristics that make it unpleasant, 
like a potential not bounded from below, we think that these characteristics 
are not so hard to be avoided. For instance, if we study the motion in the 
region around a maximum of a twice differentiable potential function then the 
approximation by a quadratic potential is valid. Particularly, this happens 
when the decay from the unstable equlibrium position is studied. 

As we have stated in the introduction, this same method has been apphed 
in the past to some models in quantum mechanics and classical statistical 
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mechanics. In the future we will try to generalize our results to more general 
models in both contexts and study the conceptual implications of the time- 
asymmetric formulation of the theory. 
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